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Abstract 

We present the results of 12 CO( J = 1-0) mapping observations toward four inter- 
acting galaxies in early and mid stages of the interaction to understand the behavior 
of molecular gas in galaxy-galaxy interaction. The observations were carried out using 
the 45-m telescope at Nobeyama Radio Observatory (NRO). We compared our CO 
total flux to those previously obtained with single-dish observations and found that 
there are no discrepancy between them. Applying a typical CO-H2 conversion factor, 
all constituent galaxies have molecular gas mass more than 10 9 M & . Comparisons to 
Hi, Ks and tracers of SF such as Ha, FUV, 8 /xm and 24 //m revealed that the dis- 
tribution of molecular gas in interacting galaxies in the early stage of the interaction 
differs from atomic gas, stars and star-forming regions. These differences are not ex- 
plained without the result of the interaction. Central concentration of molecular gas 
of interacting galaxies in the early stage of the interaction is lower than that of iso- 
lated galaxies, which suggests molecular gas is distributed off-centre and/or extends 
in the beginning of the interaction. 

Key words: galaxies: individual (Arp 84, VV 219, VV 254, the Antennae 
Galaxies) — galaxies: interactions — galaxies: ISM — ISM: molecules 



1 



1. Introduction 

Galaxy-galaxy interactions including mergers and collisions play important roles in 
galaxy formation and evolution. Close interaction with other galaxies not only changes the 
morphology and the kinematics of stars and gas in a galaxy (Toomre & Toomre 1972) but also 
enhances star-formation (SF) activity. Larson & Tinsley (1978) attribute large colour disper- 
sion of interacting systems in colour-colour diagrams compared to isolated galaxies to bursts 
of SF in interacting galaxies. Subsequent observational studies using tracers of SF such as Ha 
(Bushouse 1987; Kennicutt et al. 1987), far-infrared (FIR) (Bushouse, Werner & Lamb 1988), 
mid-infrared (MIR) (Joseph & Wright 1985; Lonsdale, Persson & Matthews 1984) and radio 
continuum (Stocke 1978; Sulentic 1976) all prove the enhancement of SF activity in interacting 
galaxies. These global and statistical studies show that a galaxy interaction makes a great 
contribution to SF activity of galaxies. 

SF occurs through a self-gravitational disruption of molecular gas clouds in interstellar 
space. It is very important to know how interstellar matter, especially, molecular gas is affected 
by interactions of galaxies to understand the mechanism of active SF in interacting galaxies. 
Since infrared (IR) luminous galaxies are expected to have strong CO emission, early CO 
observations toward extragalaxies were mainly focused on Luminous Infrared Galaxies (LIRGs) 
and Ultra Luminous Infrared Galaxies (ULIRGs) whose IR luminosity is larger than 10 11 and 
10 12 Lq, respectively (Soifer et al. 1984). Although these observations revealed that LIRGs 
have a large amount of molecular gas (> 10 10 M & ) and that there is a good correlation between 
CO luminosity and IR luminosity, the reasons for the enhanced SF activity have not been well 
understood mainly due to low spatial resolution (e.g., Young et al. 1989; Young & Scoville 
1991). 

With the improvement of receiver sensitivity, it has been possible to observe normal 
galaxies in CO and to compare SF activity between LIRGs and normal galaxies. The correlation 
between CO and IR luminosities was still found for normal galaxies and high star formation 
efficiency (SFE) was found in interacting galaxies compared to normal spiral galaxies (Young 
et al. 1986b; Sanders, Scoville & Soifer 1991; Combes et al. 1994). Considering these results, 
Solomon & Sage (1988) suggest that interacting galaxies form stars more effectively than normal 
galaxies. Recent observational studies, however, found some facts that oppose this scenario. 
Statistical studies of interacting galaxies show that molecular gas mass in interacting galaxies 
is higher than that in field galaxies (Braine & Combes 1993; Casasola, Bettoni & Galletta 
2004; Combes et al. 1994; Zhu 2001). Moreover, Casasola, Bettoni & Galletta (2004) also show 
that, contrary to the previous results, SFE is not higher in ULIRGs where star formation rate 
(SFR) is higher than field galaxies. The discrepancies of SFE in previous papers occurred 
because past observations of molecular gas were pointed at only the centre of galaxies. For 
such observations, it is difficult to get even total CO flux and a distribution of molecular gas 
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in galaxies. 

Zhu (2001) showed that the enhancement of SFR begins in the mid stage of the in- 
teraction and most active SF happens at the late stage of the interaction classified by their 
morphologies. This fact suggests that the conditions of triggering starburst have been satisfied 
in earlier stages of the interaction. The key to the onset of starburst is suggested by the numer- 
ical simulations by Barnes & Hernquist (1996). They revealed that gravitational torque makes 
interstellar gas infall into the central regions of progenitors during the interaction. Since the 
density of interstellar gas increases at the central regions due to the gas infall, SF at the galactic 
centre takes place efficiently. High resolution interferometric CO observations toward LIRGs 
clarified that molecular gas in LIRGs is concentrated in the central 1 kpc region (Scoville et al. 
1991; Yun, Scoville & Knop 1994). Since most LIRGs are regarded as being in the late stage 
of the interaction, measurements of the central concentration of molecular gas in earlier stage 
of the interaction are very important in order to explain when and how molecular gas falls into 
the nuclei. Therefore, observations of interacting galaxies in earlier stages of the interaction 
are required. 

To examine whether changes of molecular gas properties happen under the interaction, 
it is necessary to observe molecular gas which covers the entire region of interacting galaxies 
with an enough spatial resolution to resolve their main structure. The 25-BEam Array Receiver 
System (BEARS) (Sunada et al. 2000) equipped on the 45-m telescope at NRO is one of the 
most suitable instruments for such observations. Thus we conducted 12 CO(J = 1-0) mapping 
observations of interacting galaxies using the 45-m telescope. 

In this paper, we present the data from our 12 CO(J = 1-0) mapping observations toward 
four interacting galaxies in the early and the mid stages of the interaction. We investigate the 
influence of the interaction through the comparison between a central concentration of molecular 
gas and that of stars. The sample selection and the details of observations are described in 
section 2. The CO data and comparisons to previous observations are presented in section 3. 
Central concentration of molecular gas and stars are discussed in section 4. We will discuss 
molecular gas fraction and its relationship with star-forming activity in detail in forthcoming 
papers. 

2. Observations 

2.1. Target Selection 

Since our purpose is to study the properties of molecular gas in the early stages of 
the interaction, we selected interacting galaxies from the sample of Zhu (2001) based on the 
requirements as follows: 

1. Stage of the interaction 
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Fig. 1. The DSS image of our sample galaxies. The rectangles indicate our mapping area 
(3.'5 x for Arp 84, VV 219 and VV 254 and 2'.5 x 2'.5 for the Antennae Galaxies). 

To investigate the effects of the interaction on molecular gas in interacting galaxies 
and how SF is enhanced, the target galaxy pairs are in the early and the mid stages of 
the interaction. There are many classifications of interacting galaxies and most of them are 
determined by their optical morphologies. In this paper, we adopt the definition proposed by 
Zhu (2001). That is, interacting galaxies whose R < 1.5D 2 ^ were adopted as in the early stage 
of the interaction, where R is the separation between the nuclei and -D25 is the diameter of 
the primary galaxy, using the brightness contour Bt = 25 mag arcsec -2 , and showing small 
morphological disturbances. Interacting galaxies which come into contact with companions in 
projected and which have severely disturbed morphology, i.e. colliding galaxies., are classified 
as mid stage of the interaction. 

2. Intensity of CO emission 

Typically, the CO emission from the galactic disc region is a few times weaker 
than that from the central region. Therefore, we selected galaxy pairs in which CO emis- 
sion was detected from the central region with 200 Jy km/s in previous single dish observations. 

3. Size of galaxy pair 

To resolve major structures, small interacting galaxies are not suitable for our study. 
We choose galaxy pairs whose apparent size is as large as possible, for instance, 2/5x2/5. 

According to our requirements, we choose four interacting galaxies; Arp 84, VV 219, 
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VV 254 and the Antennae Galaxies as the target galaxy pairs. Basic information for each 
galaxy pair is described below. 

Arp 84 (NGC 5394 and NGC 5395) 

Arp 84 is an interacting galaxy pair in the early stage of the interaction with the 
distance of 53 Mpc and the systemic velocity of 3471 km s" 1 . Digital Sky Survey (DSS) 
image of Arp 84 is shown in Fig. 1(a). It consists of NGC 5394 (north-west) and NGC 5395 
(south-east) with a projected separation of 38 kpc and colliding in a prograde-retrograde 
interaction. Kaufman et al. (1999) suggested that NGC 5394 is nearly face-on from the narrow 
H I line width and Arp 84 has undergone a first encounter from a numerical simulation. 

VV 219 (NGC 4567 and NGC 4568) 

VV 219 galaxy pair is one of the nearest interacting galaxies located in the Virgo 
Cluster, for which we can resolve the main structures and the colliding area (see, Fig. 1(b)). 
The projected separation between two galaxies, NGC 4567 (north-west) and NGC 4568 
(south-east), of 11 kpc and their weakly disturbed optical morphology (e.g., a tidal bridge, 
tails and warped disks) implies this galaxy pair is in a very early stage of interaction, probably 
the first encounter for this interacting system. The distance of VV 219 derived from the 
systemic velocities, 2249 km s _1 , is 30 Mpc (using H = 73 km s _1 Mpc -1 ), but the distance of 
M 87 at the centre of the Virgo Cluster is given as 16 Mpc (Mei et al. 2007). This discrepancy 
may be due to the high peculiar velocity of VV 219 in the cluster potential. Since the fact of 
a projected distance between VV 219 and M 87 of 1 Mpc suggests that VV 219 is located in 
the inner region of the Virgo Cluster, we adopt 16 Mpc for the distance of VV 219. 

VV 254 (UGC 12914 and UGC 12915) 

This pair is a head-on colliding galaxy pair in which the disks of the galaxies collide 
parallel as shown in Fig. 1(c). The distance to VV 254 is 62 Mpc which is the most distant 
interacting galaxies in our targets. The galaxies located in the south-west and north-east are 
UGC 12914 and UGC 12915, respectively. The projected separation between two galaxies is 20 
kpc and the systemic velocity is 4353 km s _1 . The optical morphology of both UGC 12914 and 
UGC 12915 shows little disturbance. VV 254 is known as Taffy I and have a unique feature, 
that is, synchrotron emission is seen in the bridge between the two galaxies (Condon et al. 
1993). The direction of the linear polarization in the bridge region is aligned and perpendicular 
to the discs which suggests the trail of the collision. Since interacting galaxies VV 769 (Taffy 
II) which show the same features collide in a head-on collision (Condon, Helou & Jarrett 
2002), a head-on collision may induce aligned synchrotron emission with aligned polarization. 
The mechanism of the synchrotron emission is still not clear. Condon et al. (1993) estimate 
the time elapsed since the collision for VV 254 attributing the synchrotron emission to the 
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Table 1. Main parameters of galaxy pairs of our sample 



Name 


Velocity 


Distance 


Separation 


Resolution 




(km s ) 


/TV /T. \ 

(Mpcj 


(kpc) 


(kpc) 


I 1 ) 


\ l ) 




( i ^ 




Arp 84 


3471 


53 


38 


5.0 


VV 219 


2249 


16 


11 


1.5 


VV 254 


4353 


G2 


20 


5.8 


Antennae Galaxies 


1705 


21 


7 


2.0 


Column (1): Name oi galaxies pairs. 



Column (2): Heliocentric velocity in local standard of rest from NASA/IPAC Extragalactic Database(NED). 
Column (3): Distance of the galaxy pairs from NED. 

Column (4): Projected physical separation between two galaxies of a pair. 

Column (5): Linear spatial resolution corresponding to the effective angular resolution, 19"3, 

at the distance of galaxy pairs. 

collision and conclude that the collision occurred ~ 2 x 10 7 years ago. This means VV 254 
experienced the head-on collision recently. 

The Antennae Galaxies (NGC 4038 and NGC 4039) 

The Antennae Galaxies (Arp 244) are one of the representative prograde-prograde 
interacting galaxies and named for their elongated tidal tails. A northern galaxy is NGC 4038 
and a southern galaxy is NGC 4039. The distance to the Antennae Galaxies is 21 Mpc 
(Schweizer et al. 2008). Since they are one of the nearest interacting galaxies and show active 
SF in their distribution of molecular gas region, many observations have been carried out. 
According to the definition by Zhu, Seaquist & Kuno (2003), the Antennae Galaxies are clas- 
sified as interacting galaxies in mid stage of the interaction with a projected separation of 7 kpc. 

In table 1, basic data of our target interacting galaxies are summarized. Following the 
pair data, basic data of constituent galaxies of the target are shown in table 2. 

2.2. 12 CO(J = 1-0) Observations 

12 CO observations were made toward the target interacting galaxies with the 45-m tele- 
scope at Nobeyama Radio Observatory (NRO) 1 . We used BEARS which has twenty-five beams 
with a beam separation of 41 "1 (Sunada et al. 2000) as the receiver front-ends. The typical 
system noise temperature during the observations was T sys ~ 250-650 K (in double sideband 
mode) depending on mainly the weather. The full width at half-power (FWHP) was 16" at the 
rest frequency of 12 CO(J= 1-0)(115. 271204 GHz). 

1 NRO 45-m radio telescope is operated by Nobeyama Radio Observatory, a branch of National Astronomical 
Observatory of Japan. 
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Table 2. Main parameters of constituent galaxies of galaxy pairs 



Galaxy pair name 


Constituent galaxies name 


Right Ascension Declination 


Morphology 


Inclination Orbit 


Rk20 












(deg) 




(arcsec) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Arp 84 


NGC 5394 


13 h 58 m 33?7 


+37° 27' 13" 


SB(s)b pec 





P 


34 




NGC 5395 


13 h 58 m 37?9 


+37°24'28" 


SA(s)b pec 


58 


R 


71 


VV 219 


NGC 4567 


12 h 36 m 32f7 


+11°15'29" 


SA(rs)bc 


44 


R 


70 




NGC 4568 


12 h 36 m 34?2 


+11°14'20" 


SA(rs)bc 


58 


R 


94 


VV 254 


UGC 12914 


00 h 01 m 38?3 


+23°29'01" 


(R)S(r)cd pec 


61 


H 


55 




UGC 12915 


00 h 01 m 41f9 


+23° 29' 45" 


S? 


73 


H 


40 


Antennae Galaxies 


NGC 4038 


12 h 01 m 53?0 


-18°52'03" 


SB(s)m pec 




P 


73 




NGC 4039 


12 h 01 m 53f5 


-18°53'10" 


SA(s)m pec 




P 


81 



Name of the constituent galaxy 
Column (2): Right ascension from NED. 
Column (3): Declination from NED 
Column (4): Morphological type from NED. 

Column (5): Inclination angle. Reference: Kaufman et al. (1999) for Arp 84, this work for VV 219 and 

Giovanelli et al. (1986) for VV 254. Since the Antennae Galaxies have complex morphologies, the inclination is not fixed and is not given. 

Column (6): Orbital type judging from their morphology; P, R, H represent a prograde, retrograde and head-on collision, 

respectively. 

Column (7): The fitted radius at 20 mag arcsec -2 in the ifo-band. 

We employed the On-The-Fly (OTF) mapping technique (Sawada et al. 2008) to map the 
area of 3f5 x 3'5 which covers the whole interacting system except for the Antennae Galaxies 
for which the central 2'.5 x 2.'5 region was observed. The mapped regions for each galaxy pair 
are indicated in Fig.l with rectangles. Throughout the observations, the data were sampled 
every 0.1 s. The OFF points were taken at 7' from the centre of the interacting system. 

Twenty-five digital autocorrelators (AC45) (Sorai et al. 2000) were used as the back-ends. 
The total band width is 512 MHz (1024 channels) which corresponds to velocity coverages of 
1332 km s" 1 at 115 GHz. The frequency resolution is 500 kHz which corresponds to a velocity 
resolution of 1.3 km s _1 at 115 GHz. 

The telescope pointing accuracy was checked every hour by observing SiO maser lines 
from nearby evolved stars or 43 GHz continuum from bright QSOs with a SIS receiver (S40) 
or a HEMT receiver (H40). The antenna temperature, T\, was obtained by the chopper- wheel 
method correcting for atmospheric and antenna ohmic losses. The main-beam efficiency of the 
telescope at 115 GHz was r/ MB = 0.39 ± 0.03, 0.34 ± 0.03 and 0.31 ± 0.03 for 2004-2006, 2006- 
2007 and 2008-2009 seasons, respectively. The intensity given in this paper is the main-beam 
brightness temperature, which is defined by T MB = T^/tjmb- 
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2.3. Data Reduction 

The data were reduced with the software packages NOSTAR (Nobeyama On-The-Fly 
Software Tools for Analysis Reduction) (Sawada et al. 2008). We removed the data whose 
pointing error was larger than 5" to avoid the systematic intensity loss. Since BEARS is a 
double sideband (DSB) receiver, we calibrated the sideband ratios of each beam of BEARS 
adopting the scaling factors derived from the observations of a bright standard source using 
both BEARS and the single sideband (SSB) receiver S100. The calibration error was less than 
5 % for each array. After the calibration, we made baseline fitting to remove linear baselines. 
The data were convolved with a Gaussian-tapered Bessel function 



Ji(r/a) 

: x exp 

r/a 



r\ 2 
b) 



1.55/vr, 6 = 2.52 (1) 



to create a cube data with 7 '.'5 grid spacing and a velocity resolution of 20 km s" 1 . In the 
equation, J\ is the Bessel function of the first kind of order 1 and r the distance between data 
point and grid point in pixel. The effective angular resolution of the data is 19"3. Finally, 
we fitted and subtracted linear baselines for each spectra again. It is very important to get 
good baseline for observations of interacting galaxies because of their completely unpredicted 
distribution of molecular gas and wide velocity width of their CO line. 

3. Results 

3.1. Spectrum maps 

The spectra of 12 CO emission of Arp 84, VV 219, VV 254 and the Antennae Galaxies 
obtained with the 45 m telescope are shown in Fig.2 — 5. 

Achieved sensitivities were calculated by averaging the rms noise per channel of the 
spectra in emission free regions and are summarized in Table 3. The uncertainties in the 
integrated intensity were estimated using the next equation; 



AI = ^AP n + AP h (2) 

where 

AI n = T vms dv.i(dv c /dvi) 1/2 (3) 

AI b = T rms dv t (dv c /dv b ) 1/2 (4) 

and AI n the rms error in the spectrum, AI b the error on the determination of the baseline. 
In these equations, T rms means the rms noise per channel and dvi is velocity range over which 
the spectrum was integrated, dv c the smoothed channel width and dv b the velocity width over 
which the baseline was fitted. Typically, AI n is less than 10 % of the total intensity and is 
mainly weather dependent. The baseline uncertainty AI b is usually 5-10 %. 
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Fig. 2. Spectrum map of CO(1-0) emission of Arp 84. The horizontal axis and verti- 
cal axis for each spectra are velocity in local standard of rest and main beam tempera- 
ture. Mapped region is 210" x 210" which corresponds 54 kpc x 54 kpc in linear scale. 
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Table 3. Achieved 12 CO(1-0) sensitivities 



Source 


rms 


Observing time^ 




(mK) 


(hours) 


Arp 84 


14 


32 


VV 219 


19 


40 


VV 254 


25 


24 


Antennae Galaxies 


37 


10 



* rms noise per channel with 20 km s 1 velocity resolution, 
t the net observing time after subtracting bad data. 
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Fig. 4. Spectrum map of CO(l-O) emission of VV 254. Mapped re- 
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Fig. 5. Spectrum map of CO(l-O) emission of the Antennae Galaxies. Mapped 
region is 150" x 150" which corresponds 15 kpc x 15 kpc in linear scale. 



3.2. Comparisons with previous single-dish observations 

Past CO observations toward our sample were used to check our intensity calibration. 
At first, we derived the flux from our 12 CO integrated intensity map. The formula to derive 
the flux density (Sqo m Jy km -1 ) per pixel in maps from the brightness temperature (2b i n 
K) is as follows: 

2k 



CO 



A 2 



(5) 



where k is the Boltzmann constant, A is the wavelength of the observed 12 CO(1-0) line in rest 
frame (2.6 mm) and Q is the solid angle of a pixel of our map 2 . We adopted a conversion factor 
Jy K -1 pixel -1 is 0.5 and we use T MB instead of T B . 

The total CO flux of Arp 84 derived from our map including calibration uncertainties is 
243 ± 38 Jy km s -1 for NGC 5394 and 406 ± 63 Jy km s" 1 for NGC 5395. These are consistent 



2 The solid angle occupied by a pixel (7"5 x 7"5) is 1.32 x 10 1 steradian 
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with the results of Zhu et al. (1999) obtained with the NRAO 12-m telescope (201.3 ± 16.7 
Jy km s" 1 for NGC 5394 and 547.8 ± 55.2 Jy km s" 1 for NGC 5395) within the error of their 
flux estimation (20-30 %). Zhu et al. (1999) derived the total flux of NGC 5395 by assuming 
exponential radial distribution of molecular gas. The real radial distribution appears steeper 
than the exponential model as seen in our map (Fig. 6 (a)). This may explain why the total 
flux of NGC 5395 by Zhu et al. (1999) is larger than ours. Kaufman et al. (1999) derived the 
total flux of 177 Jy km s" 1 for NGC 5394 and 203 Jy km s" 1 for NGC 5395 using the 20-m 
telescope at Onsala Space Observatory. Since they observed only 13 positions and they are not 
performed with Nyquist sampling, the flux is lower than our results. 

The total flux of VV 219 system derived from our map is 714 ± 111 Jy km s _1 for 
NGC 4567 and 1664 ± 258 Jy km s" 1 for NGC 4568. Total CO flux of NGC 4567 and 
NGC 4568 reported by Young et al. (1995) are 500 Jy km s _1 and 1050 Jy km s _1 , respectively. 
Since they observed only one position and their beam does not cover the entire region of the 
galaxies, lower flux is expected. Although Chung et al. (2009) mapped the entire region of the 
galaxy pair using the Five College Radio Astronomy Observatory (FCRAO) 14-m telescope, 
the total flux they derived is 50 % lower than ours (1500 ± 170 Jy km s _1 ). Since they found 
their flux is ~25 % smaller than that obtained by BIMA (Heifer et al. 2003) and attributed 
the discrepancy to the systematic error in calibration, their values may be underestimated. 

The total CO flux of VV 254 derived from our map is 1422 ± 220 Jy km s" 1 . This 
is consistent with the results obtained with the IRAM 30-m telescope by Braine et al. (2003) 
(1246 Jy km s -1 ) and those with the NRAO 12-m telescope by Zhu, Seaquist & Kuno (2003) 
(1012.1 ± 171.6 Jy km s -1 ), taking into account that our mapped region is larger than theirs. 

Our total CO flux for the Antennae Galaxies is 3153 ± 426 Jy km s -1 , which is 10 % 
lower than that derived by Gao et al. (2001) (3595 Jy km s _1 ) using the NRAO 12-m telescope, 
consistent with the FCRAO 14-m telescope (Young et al. 1995) (2844 Jy km s" 1 ), and 30 % 
higher than that detected by Zhu, Seaquist & Kuno (2003) (2370 Jy km s _1 ). The observed 
region of Zhu, Seaquist & Kuno (2003) is smaller than our map although both observations were 
accomplished by the NRO 45-m telescope. We confirmed that our spectra are consistent with 
those in Zhu, Seaquist & Kuno (2003) within their observed region. On the other hand, Gao 
et al. (2001) observed a 4' x 4' region which is larger than ours (2'.5 x 2'5) and weak emission 
is seen in the region where we did not observe. Therefore, the main cause of the difference of 
the total flux in these observations is explained by the difference of the mapping areas. 

We derived the molecular gas mass of each galaxy pair from CO flux using the following 
formula (Kenney & Young 1988): 

M H2 = 3.9 x 10~ 17 Acq D 2 S co (6) 

where M Ha is molecular gas mass in the unit of M , X co the N H2 //co conversion factor and 
D distance to the galaxy in Mpc. We adopt a standard conversion factor X C o = 1-8 x 10 20 
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cm -2 (K km s -1 ) -1 (Dame, Hartmann & Thaddeus 2001). On deriving molecular gas mass, we 
divide the galaxy pair into two constituent galaxies in the manner described below. We fit two 
ellipses to reproduce the contour of 20 mag arcsec -2 in ifs-band. We regard the region inside 
of fitted ellipses as the extent of the galaxy. For VV 219 and the Antennae Galaxies, since 
fitted ellipses overlap, we divide two galaxies with a line which crosses the points where two 
ellipses intersect. The bridge region of VV 254 is defined as Braine et al. (2003). Fitted ellipses 
are drawn in red in Fig. 6(d), Fig. 7(d), Fig.8(d) and Fig. 9(d). Total CO flux and molecular gas 
mass of each galaxy are summarized in table 4. All galaxies in our sample have molecular gas 
mass of more than 10 9 M . 

It needs to notice, however, that Aco varies between galaxies (Nakai & Kuno 1995) and 
is not constant even within a galaxy (Casasola et al. 2007). Observations of nearby galaxies 
revealed that Aco is a function of metallicity (Arimoto, Sofue & Tsujimoto 1996; Wilson 
1995; Leroy et al. 2011). Most important fact of Aco for our studies that (U)LIRGs most 
of which are interacting galaxies in the late stage of the interaction have lower Aco than our 
Galaxy (Solomon et al. 1997; Downes & Solomon 1998). Although past studies attribute low 
conversion factor to warmer molecular gas temperature and larger CO line width and velocity 
dispersion (Maloney & Black 1988; Downes & Solomon 1998; Narayanan et al. 2011), the 
origin of low Aco in interacting galaxies is still unknown. Aco factor in the Antennae Galaxies 
and VV 254 are also derived. The Antennae Galaxies have an one order lower Aco than the 
standard conversion factor (Zhu, Seaquist & Kuno 2003). Zhu et al. (2007) found that A C o in 
VV 254 is about 4 times lower compared to the standard conversion factor. This fact suggests 
that estimated molecular gas mass may be overestimated even for Arp 84 and VV 219 for which 
we do not have available Acq- 
s'. 3. Distributions of molecular gas and comparison with other Wavelength 

We present distribution of molecular gas and velocity field of our sample galaxy pairs 
and compare between the distribution of molecular gas and those of H I, ifs-band, Ha or FUV 
and 8 /im or 24 /jm. H I data which are obtained with VLA by Iono, Yun & Ho (2005) except for 
the Antennae Galaxies which was taken by Hibbard et al. (2001) and i^s-band 2MASS archival 
data (Jarrett et al. 2000) are used for the comparison. Ha data of VV 219 is acquired from 
Koopmann, Kenney & Young (2001) and for the Antennae Galaxies from Xu et al. (2000). FUV 
data are taken from GALEX data archive, and 8 /im and 24 /im are from Spizter IRAC and 
MIPS data archives (Smith et al. 2007). The angular resolution of Hi is summarized in Table 
5. Point spread functions of Ks, Ha, FUV, 8 /im and 24 /im are typically 2", 2", 4", 2" and 
5", respectively. For each galaxy pair, integrated intensity map, velocity field, Hi integrated 
intensity map and ife-band image, Ha, FUV, 8 /im and 24 /im are shown in Fig. 6 — 9. 
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Table 4. Derived Properties of Molecular Gas 



Name 


Total Flux (Jy km s^ 1 ) 


M H2 (1O 9 M ) 


c 20 t 


Cj 


Arp 84 


649 ± 103 


12.76 ± 2.00 






NGC 5394 


243 ± 38 


4.77 ± 0.75 


0.42 


0.80 


NGC 5395 


406 ± 63 


7.99 ± 1.25 


0.28 


0.43 


VV 219 


2378 ± 369 


4.26 ± 0.66 






NGC 4567 


714 ± 111 


1.28 ± 0.20 


0.49 


0.46 


NGC 4568 


1664 ± 258 


2.99 ± 0.46 


0.48 


0.62 


VV 254 


1422 ± 220 


38.4 ± 5.94 






UGC 12914 


380 ± 59 


10.3 ± 1.60 


0.36 


0.68 


UGC 12915 


554 ± 86 


14.9 ± 2.31 


0.36 


0.79 


the Bridge 


488 ± 75 


13.2 ± 2.03 






Antennae Galaxies 


3153 ± 426 


9.33 ± 1.26 






NGC 4038 


2080 ± 281 


6.15 ± 0.83 


0.45 


0.48 


NGC 4039 


1073 ± 145 


3.18 ± 0.43 


0.46 


0.50 



t Central concentration of molecular gas C20 is defined by equation 7. 
t Central concentration of star C* is defined by equation 8. 



Table 5. Angular Resolution of Hi data used for comparison with our 12 CO(1-0) observations 

Galaxy Pair Angular Resolution (arcsec) 
Arp 84 17.8 x 16.3 

VV 219 19.9 x 14.0 

VV 254 18.0 x 18.0 

Antennae Galaxies 11.4 x 7.4 



3.3.1. Arp 84 

Fig. 6 (a) shows the CO integrated intensity map of Arp 84. The CO map shows that 
NGC 5394 has a strong peak near the centre. There is an offset between the position of the 
peaks of CO and Ks in NGC 5394. The CO peak is located about 5" west of the Ks peak. 
Although the offset is smaller than our beam size, the offset is also seen in the map obtained 
with an interferometer (Iono, Yun & Ho 2005). The distribution of CO is slightly elongated 
from east to west in NGC 5394, while the Ks image shows the tidal tail from north to south. 
Hi is very weak in NGC 5394, while a strong peak is seen in CO. This suggests that most of the 
interstellar gas exist as molecules in NGC 5394. On the other hand, Fig. 6 shows that NGC 5395 
has extended structure in all data. Ks image shows a strong peak at the centre and a long tidal 
arm from south to north. H I and CO is distributed along the tidal arm. H I also has weak tidal 
arm at the northern part of NGC 5395 where no CO emission is seen (Declination > 37°27'00"). 
There are H I and CO peaks at the northern end of the tidal arm, while there is no peak of CO 
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Fig. 6. The crosses represent the galactic centre of NGC 5394 and NGC 5395. The contours are 2.81 
x 2, 3, ... K km s" 1 for all images except for the image (d) where contours are 3160 + 40 x 1, 2, 
3, ... km s _1 . (a) 12 CO(1-0) integrated intensity map and contours, (b) 12 CO(1-0) velocity field, (c) 
12 CO(1-0) contours superposed on the Hi image (Iono, Yun & Ho 2005). (d) 12 CO(1-0) contours su- 
perposed on the ifs-band image (Jarrett et al. 2000). Red eclipses, result of the fit whose method 
is described in chapter 3.2. (e) 12 CO(1-0) contours superposed on the GALEX-FUV image (Bianchi 
et al. 2003). (f) 12 CO(1-0) contours superposed on the MIPS 8 micron image (Smith et al. 2007). 



and H I at the centre of NGC 5395 unlike NGC 5394. Both NGC 5394 and NGC 5395 show the 
discrepancy between CO /H I and Ks. The discrepancy may imply the different behavior of gas 
and stars in the interaction. While the velocity field of NGC 5394 indicates that the galaxy is 
almost face on, NGC 5395 shows rotating features and weakly disturbed at the northern end 
of the tidal arm. 
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Fig. 6(e) and (f) show the distribution of FUV and 8 /im both of which trace SF activity. 
In NGC 5394, the distribution of FUV shows that SF activity is very high at the centre of the 
galaxy. Moreover, star-forming regions are also found along the tidal spiral arms seen in Ks. 
There is a local peak of SF in the middle of northern tidal arm of NGC 5394. On the other 
hand, NGC 5395 shows an extremely complex distribution of star-forming regions. Although 
a peak of SF activity is located near the peak of CO, it does not coincide with the CO peak. 
The distribution of star-forming regions extends toward NGC 5394 and is distributed all over 
NGC 5395. The distribution of star-forming regions traced by FUV shows that the tidal arm 
of NGC 5395 is connected with the northern tidal tail which is not seen in the distribution of 
CO and ife-band. The southern tidal arm of NGC 5395 traced by FUV extends to NGC 5394 
and has its strongest peak in this galaxy. 

The distribution of 8 /im of NGC 5394 shows a peak at the centre and in tidal spiral 
arms as seen in the FUV map. Weak 8 /im emission from tidal spirals of NGC 5394 compared 
to FUV implies that the dust extinction at the tidal spirals of NGC 5394 is small. The fact 
that star-forming regions traced by both FUV and 8 /im are seen at the centre of NGC 5394 
suggests strong SF occurs and is concentrated to the centre of NGC 5394. The distribution of 
8 /im of NGC 5395 resembles to i^s-band. Each local peaks of 8 /im in NGC 5395 are located 
near the CO peaks. 

Discrepancies of distribution between interstellar gas (i.e., CO and Hi) and star-forming 
regions (i.e., FUV and 8 /im) are also seen in field spiral galaxies, such as M 100 (Sempere 
& Garcia-Burillo 1997), NGC 5194 (Calzetti et al. 2005), IC 342 (Hirota et al. 2010) and 
NGC 3147 (Casasola et al. 2008). In field spiral galaxies, the offsets between ISM and star- 
forming regions may be explained with timescale variations and SFE in the result of a spiral 
density wave. Molecular gas forms stars in the density wave and then newly formed stars are 
separated from their parent molecular cloud. While molecular gas moves to other condensation 
along the density wave, the young stars emit FUV and can heat only nearby dust. In this 
scenario, CO peaks are regarded as the future SF site and decouple from star-forming tracers. 
The inconsistency of distribution between ISM and star-forming regions in NGC 5395 is more 
eminent than that in spiral galaxies which is previously reported. This fact may be explained 
that the disturbance and accumulation of ISM due to the interaction in NGC 5395 is much faster 
than SF timescale and/or SF timescale becomes longer than the field spiral galaxies. However, 
latter possibility is not likely to happen, since previous observations found that SFE, which is 
an inverse of SF timescale, in interacting galaxies is higher or at least the same compared to 
that in field spiral galaxies (Young et al. 1989; Casasola, Bettoni & Galletta 2004). In either 
case, the CO peaks found in the edge of the northern tidal tale in NGC 5395 would be the 
birth site of new stars. 
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3.3.2. VV 219 

Fig. 7(a) shows the CO integrated intensity map of VV 219. The figure shows that both 
galaxies have a central peak as often seen in nearby spiral galaxies (e.g., Kuno et al. 2007). 
Although the distribution of molecular gas in these galaxies is not significantly disturbed by 
the interaction, the distribution of molecular gas in NGC 4568 shows asymmetry with respect 
to the minor axis of the galaxy. That is, the extension of distribution of molecular gas on the 
northern side is larger connecting to NGC 4567 than the southern side. The molecular gas 
smoothly connects from the northern edge of NGC 4568 to the eastern edge of NGC 4567. 
The i^s-band image (Fig. 7(d)) shows spiral structure of both galaxies (Jarrett et al. 2003). 
Although the angular resolution of CO is not enough to resolve the spiral arms, we can see 
that the molecular gas tends to concentrate on the spiral arms. NGC 4567 is known as a Hi 
deficient galaxy (Koopmann & Kenney 2004). In Fig. 7(c), most of the detected Hi emission 
in NGC 4567 is distributed within the lowest contour of the CO map which corresponds to 
the surface density of molecular gas of 19 M Q pc -2 . CO and Hi extend toward the opposite 
direction in NGC 4568. That is, while CO extends toward the northern side, H I extends toward 
the southern side. The strongest Hi peak is located on the elongated CO structure toward the 
north. 

The velocity field of molecular gas is shown in Fig. 7(b). It shows kinematics of molec- 
ular gas more clearly as compared with the previous studies whose field of view or angular 
resolution were limited (e.g., Iono, Yun & Ho 2005). The velocity fields of both galaxies are not 
significantly disturbed largely and overlap smoothly. The radial velocity increases from west to 
east for NGC 4567 and from south to north for NGC 4568. If we assume that NGC 4567/68 
possess trailing arms, both galaxies rotate clockwise on the celestial sphere and are undergoing 
a retrograde-retrograde collision. Undisturbed but smooth connection at the distribution of 
molecular gas region in the velocity field confirms that NGC 4567 and NGC 4568 are not an 
apparent collision but are really colliding. 

As shown in Fig. 7(e) and (f), there are some similarities between Ha and 24 /xm distri- 
butions. Massive star-forming regions are concentrated in the central region and spiral arms 
in both galaxies. Although most of the star-forming regions are found along the spiral arms 
seen in Ks image (R < 45"), there are some large star- forming regions on the northern side 
of NGC 4568 which corresponds to the overlap region with NGC 4567, where strong Hi peak 
and the CO elongation are seen. These results imply that the interaction of NGC 4567 and 
NGC 4568 affected the SF in the overlap region. 

3.3.3. VV 254 

The CO integrated intensity map of VV 254 is shown in Fig. 8 (a). Molecular gas in 
VV 254 is stretched out to the bridge region which is between the galaxies and its peaks 
are lopsided to the bridge region made by the head-on collision. Although the apparent size of 
UGC 12915 is smaller than UGC 12914, a large amount of molecular gas belongs to UGC 12915. 
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Fig. 7. Same as Fig. 6 but for VV 219 and (e) is 12 CO(1-0) contours superposed on the Ha im- 
age and (f) is 12 CO(1-0) contours superposed on the MIPS 24 micron image. The contours are 
7.20 x 1, 2, 3, ... K km s _1 for all images except for (b) where the contours are 2100 + 20 
x 1, 2, 3, ... km s~ 1 . The crosses represent the galactic centre of NGC 4567 and NGC 4568. 



The CO integrated intensity map also reveals that molecular gas has concentrated on the eastern 
side of UGC 12915. 

Fig.8(b) shows that kinematics of VV 254 is severely disturbed, especially in the bridge 
region. Although both UGC 12914 and UGC 12915 are nearly edge-on, the velocity field does 
not show rotation. 

Fig. 8(c) shows that the distribution of atomic gas of VV 254 is much different from that 
of molecular gas. It extends from north to south of VV 254 and most of the atomic gas is found 
in the middle of the bridge region. Atomic gas is more rich in UGC 12914 than in UGC 12915 
whereas molecular gas is more abundant in UGC 12915. Although the central region of two 
galaxies are not disturbed in Ks, the edges of the discs of the galaxies are warped toward the 
counterpart galaxy. 

Fig. 8(e) and (f) represents most star- forming regions in VV 254 are embedded in the 
discs of the galaxies and there is no star-forming region in the bridge region except for around 
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Fig. 8. Same as Fig.6 but for VV 254 and (e) is 12 CO(1-0) contours superposed on the GALEX-FUV 
image and (f) is 12 CO(1-0) contours superposed on the MIPS 24 micron image. The contours are 6.63 x 
2, 3, 4, 5... K km s _1 for all images except for (b) where the contours are 4100 + 20 x 1, 2, 3, ... km 
s _1 . The crosses represent the galactic centre for UGC 12914 (bottom right) and UGC 12915 (top left). 

the CO peak of UGC 12915. In UGC 12914, the strongest SF activity traced by FUV shown 
in Fig. 8(e) is found at the northwest disc where CO emission is weak. Similarly, the most 
active star-forming region of UGC 12915 is located about 10" from the centre and it does not 
correspond to the CO peak. Although both galaxies show weak SF activity in the centre of the 
galaxy in the FUV image, 24 /im image (Fig. 8(f)) reveals that large star- forming regions are 
obscured in the centre of both galaxies. UGC 12914 has three massive star-forming regions: 
the centre of the galaxy and both edge of the disc. In the same way as FUV, the strongest 
star-forming region traced by 24 /im is seen at the western part of the disc of UGC 12915. 

While molecular gas is plenty (a mass of more than 10 10 M@) in the overlap region 
where large amount of atomic gas is found, there is no star-forming region there. Since VV 254 
is a head-on collision galaxy and is only 10 7 years after first encounter, almost all molecular 
gas could be gravitationally unbound. The fact that the overlap region has higher 12 CO (2- 
1)/ 12 CO(1-0) and 12 CO(3-2)/ 12 CO(l-0) ratio than UGC 12914 and UGC 12915 (Zhu et al. 
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Fig. 9. Same as Fig. 6 but for the Antennae Galaxies and (e) is 12 CO(1-0) contours superposed on the 
Ha image and (f) is 12 CO(1-0) contours superposed on the MIPS 24 micron image. The contours arc 
23.7 x 1, 2, 3, ... K km s _1 for all images except for (b) where the contours arc 1500 + 20 x 1, 2, 
3, ... km s _1 . The crosses represent the galactic centre for NGC 4038 (top) and NGC 4039 (bottom). 



2007) suggests that denser and warmer molecular gas which is more direct future SF tracer is 
embedded to the overlap region. Therefore, VV 254 will produce stars at the overlap region in 
the future. 

3.3.4- The Antennae Galaxies 

Although a large amount of molecular gas exists in the centre of both galaxies, the 
largest concentration of molecular gas is seen in the overlap region (Fig. 9 (a)). The mass of 
the concentration is more than 10 9 M & . We discovered molecular gas at the north east region 
around (a, 8) = (12 h l m 57 s , -18°52'10"), which had not been found in the previous observations 
and belongs to the overlap region judging from Fig. 5, although the amount of molecular gas in 
these structures is small. On the other hand, atomic gas traces the tail structures and spiral 
arms of NGC 4038. Hi cavities are found at the centre of NGC 4038 and the overlap region 
where most molecular gas is found. Additionally, the peaks of H I gas is located in at the north 
of the overlap region, the south of NGC 4039 and the tidal spiral arm of NGC 4038. These 
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regions do not correspond to the local peaks of CO emission. Fig. 9(d) shows that old stars have 
smooth connections from NGC 4038 to NGC 4039. Most molecular gas is distributed within 
the radius where brightness of Ks-b&nd is 20 mag arcsec -2 in the discs of these galaxies. The 
velocity field shows complex kinematics of molecular gas (Fig.9(b)). While NGC 4039 shows 
nearly rotating features, the velocity field of NGC 4038 does not show clear kinematics. 

Ha image in Fig. 9(e) shows that there are many active star- forming regions whose lumi- 
nosity is comparable to that in the central region of both galaxies. These star-forming regions 
are distributed along the structures seen in Ks-b&nd (Fig.9(d)) and SF occurs not only in the 
region where CO emission is strong but also where CO is not strong like the northern distri- 
bution of molecular gas region. Most star-forming regions traced by 24 /im is also found in 
Ha image, although many star-forming regions are found in Ha. Besides the distribution of 
star-forming regions, their relative luminosity differs between Ha and 24 fim. Although the 
distribution of star-forming regions traced by 24 /im coincides with CO distribution, the most 
active star-forming region is found in the southern distribution of molecular gas region which 
is located about 20" (2 kpc in a linear scale) southeast from the CO peak. Taking into account 
for a discrepancy of distributions between CO and Hi like NGC 5395, the birth site of stars is 
expected to move to the northern part of the distribution of molecular gas region. 

4. Discussion 

In this section, we discuss central concentration of molecular gas in interacting galaxies 
by comparing them to isolated spiral galaxies to investigate the nature of the distribution of 
molecular gas under the interaction. 

For the comparisons, the control sample of isolated spiral galaxies was obtained from 
Nobeyama CO atlas (Kuno et al. 2007). This atlas is the CO mapping survey of 40 nearby 
spiral galaxies using the NRO 45-m telescope. Since the atlas was obtained with the same 
telescope as our data, the atlas data can be used for a direct comparison with our data. To 
remove the environmental effect, we excluded the galaxies which belong to the Virgo Cluster 
and the Coma Cluster from the atlas sample. We also removed galaxies whose Hi data is 
not available. The final control sample selected with these requirements comprises 20 galaxies 
summarized in Table 6. 

One of the scenarios which explains the explosive SF activity found in interacting galaxies 
in the late stage of the interaction is gas infall toward the galactic centre as shown by Barnes 
& Hernquist (1996). The gas infall leads to an enhancement of the density of molecular gas 
in the central region of the galaxies and active SF will be triggered using the dense molecular 
gas. In fact, high molecular gas concentration has been found in late stage interacting galaxies 
such as Arp 220 and NGC 520 whose SF is extremely active (Scoville et al. 1986; Sanders et al. 
1988). 

On the other hand, for early and mid stage interacting galaxies, Iono, Yun & Ho (2005) 
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Table 6. Control sample of isolated galaxies 



Name 


Morphology 


Velocity 


Distance 


Inclination 




D t 
WK20 


S~1 


S~1 

L * 






(km s -1 ) 


/1\ /T„^,\ 

^ivipc J 


(deg) 


(10 M J 


(arcsec) 






I 1 ) 


( 2 ) 


(3) 


/ 1 \ 

(4) 


(5) 


(6) 


{<) 


(8) 


(9) 


MGC 253 


OAT)/ \ 

bAB(s)c 


227 


2..") 


75 


1.0 


630.2 


0.78 


0.77 


Mallei 2 


HAD/ „ 

bAb(rs)bc 


-24 


O O 

2.8 


i 2 


A A 

0.9 


286.4 


A AO 

0.48 


a a a 

0.64 


IC 342 


bAB(rs)cd 


28 


3.3 


31 


2.7 


372.0 


0.47 


0.59 


UGC 2855 


bABc 


1207 


10.4 


63 


20.3 


114.3 


0.53 


0.73 


MGC 2903 


A T"> / \ "U 

bAB(rs)bc 


549 


6.3 


67 


1.6 


163.0 


0.54 


0.62 


MGC 3184 


O A T~> / \ J 

SAB(rs)cd 


594 


8.7 


21 


1.1 


114.6 


0.38 


0.57 


MGC 3351 


Or) / M 

bB(r)b 


778 


10.1 


40 


1.1 


116.3 


0.45 


0.85 


INGC 35zl 


OAT - )/ ~ \ U ~ 

bAB(rs)bc 


tao 

792 


T o 

/ .2 


/ " > 

Do 


O A 

2.9 


1 C A A 

lo4.4 


A f O 


a or* 

0.86 


INLrC OOZY 


bAB^sjb 


71 r 

715 


ii i 
11.1 


52 


o o 
O.O 


1 O C A 

loo. (J 


A £T 

O.o7 


A 7C 
0. (0 


AT/"*/"* oroi 

MGC 3631 


O A /^.\ ~ 

bA(s)c 


1154 


1 o 

21. 


1 i 


3.0 


OA A 

80.9 


A CO 

0.52 


A T£? 

0.76 


IN GO 4051 


HAD/ „\1 

bAB(rs)bc 


725 


1 T A 

17.0 


1 A 

49 


O T 

2.7 


1 AO £? 

102.6 


A A A 

U.44 


A TCT 

0.75 


INGC 410z 


O A TO /.- \"U 

bAB(s)b 


853 


1 T A 

17.0 


56 


O O 

2.2 


08.6 


A TO 

0.72 


A OA 

0.80 


AT/^1 O /IT Of? 

MGC 473b 


/T3\0 A / 

(KJbA(r)ab 


O 1 T 

317 


4.3 


1 A 

40 


A 1 

0.4 


1 TO O 

172.3 


A A 1 

0.91 


A OT 

0.87 


at/"*/"* irAirc 

MGG 5055 


bA(rs)bc 


AO 

503 


T O 

7.2 


01 


O A 

2.9 


OA A O 

204.2 


A AO 

0.48 


A TIT 

0.75 


t\t/^1 POO/"* 

MGG 5236 


O A T~> /~ \ „ 

bAB(s)c 


514 


■I r 

4.5 


24 


O A 

2.0 


312.4 


A O T 

0.87 


A TA 

0.70 


NGC 5248 


SAB(rs)cd 


1165 


22.7 


40 


8.8 


111.8 


0.68 


0.70 


NGC 5457 


SAB(rs)cd 


255 


7.2 


18 


3.2 


236.3 


0.44 


0.80 


NGC 6217 


(R)bA(rs)bc 


1355 


23.9 


34 


1.8 


73.2 


0.69 


0.57 


NGC 6946 


SAB(rs)cd 


60 


5.5 


40 


4.0 


252.5 


0.57 


0.46 


NGC 6951 


bAB(rs)bc 


1425 


24.1 


30 


5.9 


115.6 


0.59 


0.62 



* Velocity in local standard of rest (Kuno ct al. 2007). 
t Molecular gas mass from Kuno et al. 2007. 

i The radius at 20 mag arcsec -2 in the ife-band (Jarrett et al. 2003). 
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claimed that molecular gas is not concentrated to the nucleus of progenitors as expected in 
numerical simulations. Their observations, however, were carried out with the interferometer 
and mapped only central regions of the sample galaxies. The missing flux problem of the inter- 
ferometer disables measurements of the degree of concentration of molecular gas. In addition, 
because of the limited observing area, it is difficult to know the actual distribution of molecular 
gas in the early and the mid stage of the interaction in which molecular gas is distributed 
over a relatively wide area in the galactic discs. Hence, mapping observations of the whole 
region of the galaxy pair with a single-dish telescope are important to quantify concentration 
of molecular gas. 

In this paper, molecular gas concentration is defined by 
M K2 (R K2Q /2) 

° 2 ° ~ M H2 (R K20 ) ' (7) 
where M^ 2 {R) is molecular gas mass within a radius R and Rk2$ is the radius at 20 mag 
arcsec -2 in Ks-band. We took Rk2o value from 2MASS Large Galaxy Atlas (Jarrett et al. 
2003). For the galaxies which are not included 2MASS Large Galaxy Atlas, we derived Rk2o 
from the Two Micron All Sky Survey (2MASS) Ks-band archival data for isolated spiral galaxies 
(Jarrett et al. 2000) 3 . Rk2o of interacting galaxies are taken from the major axis of the ellipse 
which is defined in §3.2. For VV 219 and the Antennae Galaxies, since 20 mag arcsec -2 contour 
connects between two constituent galaxies, we exclude the regions outside the red line which 
divide two galaxies to calculate molecular gas concentration. The radius at 20 mag arcsec -2 in 
i^s-band for each interacting galaxy is shown in Fig.6(d) — Fig. 9(d) with white ellipses. C 2 q is 
listed in Table 4 and Table 6. 

Fig. 10 (a) is the histogram of C 2 q for the interacting galaxies and the control galaxies sam- 
ple. The median of C 2 q of interacting galaxies is 0.44 and the maximum is 0.49, while the median 
of the control sample is 0.56 and the maximum is 0.91. We performed the Kolmogorov-Smirnov 
test between two samples and found that the probability that both samples are extracted from 
the same parent sample is 0.02. As mentioned in §3.3, C 2 q of NGC 5394 may be underestimated 
because CO extension is relatively small as compared to other galaxies. Therefore we performed 
the Kolmogorov-Smirnov test assuming the C 2 q value of 0.9 for NGC 5394 which is comparable 
with the maximum value of the control sample. Even in this case, the possibility is 0.08. This 
result implies that molecular gas is not accumulated toward the centre of interacting galaxies 
directly by the interaction and that there is a phase in which distribution of molecular gas is 
significantly disturbed. 

For comparison, we also derived central concentration of stars, C*, using a Ks-band 

3 This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project 
of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute 
of Technology, funded by the National Aeronautics and Space Administration and the National Science 
Foundation. 
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Fig. 10. (a) Histogram of C20 for the interacting galaxies sample (upper panel) and the control sample 
(lower panel), (b) Histogram of C* for the interacting galaxies sample (upper panel) and the control 
sample (lower panel). 

image which traces stellar mass. In the same manner as C20, is defined as 
_ L*(R K20 /2) 

* ~ t (n — T' W 

where L* is Ks-b&nd luminosity. is listed in Table 4 and Table 6. The median of C* of 
interacting galaxies is 0.56 and the maximum is 0.79, while the median of the control sample is 
0.74 and the maximum is 0.87. The result of the Kolmogorov-Smirnov test of 0.20 shows these 
two samples are from same parent, although the control sample has higher value of C*. 

These comparisons of C20 and C* imply that both molecular gas and stars are disturbed 
by the interaction and that the interaction has a stronger effects on molecular gas than stars in 
the early stage of the interaction. Our observational finding that distributions of molecular gas 
in the early stage of the interaction are larger than isolated spiral galaxies is consistent with 
recent simulations which treat feedback from SF and precise physics of interstellar medium 
(Bournaud et al. 2010; Saitoh et al. 2009; Teyssier, Chapon & Bournaud 2010), while previous 
simulations which adopt isothermal interstellar medium (Barnes & Hernquist 1996; Mihos & 
Hernquist 1996) do not show a significantly disturbed distribution of gas. This fact indicates 
that direct comparisons can be done between observation and simulation about not merely of 
distribution of gas but SF associated with the change of physical properties which are discussed 
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in upcoming papers. 



5. Summary 

We made CO mapping observations of interacting galaxies in the early and the mid 
stages of the interaction with the 45-m telescope at NRO. Our high quality CO maps reveal that 
distribution of molecular gas differs from atomic gas, old stars and SF tracers suggesting the 
strong influence of the interaction even in the early stage of the interaction. Some interacting 
galaxies show that molecular gas is lopsided to the area between two constituent galaxies. 
Furthermore, we compared central concentration of molecular gas between interacting galaxies 
and isolated galaxies and found that the central concentration of molecular gas in interacting 
galaxies is lower than that in isolated galaxies. Old stars also show interacting galaxies have 
lower central concentration than isolated galaxies but the difference is not so clear as compared 
with that of molecular gas. The results suggest that even in the early stage of the interaction the 
distribution of molecular gas is strongly disturbed by the interaction and there is a phase when 
the central concentration of molecular gas decreases before gas infall toward the galactic centre. 

We would like to thank the all the staff of the NRO for observational support. This 
research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated 
by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the 
National Aeronautics and Space Administration. YT is supported by JSPS Grant-in-Aid for 
Research Activity Start-up (no. 23840007). 
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